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Abstract 

We present a preliminary measurement of CP-violating parameters S and C from fits of the 
time-dependence of meson decays to rj'K^. The data were recorded with the BABAR detector at 
PEP-II and correspond to 227 x 10^ BB pairs produced in e^e annihilation through the T{AS) 
resonance. From a maximum likelihood fit we measure the CP-violation parameters S = 0.27ib0.14 
(stat) ± 0.03 (syst) and C = -0.21 ± 0.10 (stat) ± 0.03 (syst). 
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1 Introduction 



Measurements of time-dependent CP asymmetries in meson decays through a Cabibbo-Kobayashi- 
Maskawa (CKM) favored h ccs amphtude ^ have provided a crucial test of the CKM mechanism 
of CP violation in the Standard Model (SM) [2]. Such decays to a charmonium state plus a 
meson are dominated by a single weak phase. Decays of mesons to charmless hadronic final 
states, such as (/>K°, rj'K^, K^K~K^^ K^tt^ and /o(980)-ftr'^, are expected to be dominated by pen- 
guin diagrams. If we neglect CKM-suppressed amplitudes, these decay modes have the same weak 
phase as the charmonium decays in the SM. Thus the time-dependent asymmetry measurement 
for these decays should yield an alternative measurement of sin 2/3 0. 

The processes shown in Fig.^b)-(d) are relevant for the decay B^ r]' , and there are similar 
diagrams for B~^ r]'K~^. 

All of the amplitudes for these processes have CKM suppression, but the tree diagram for B^ 
shown in Fig. D^b) is expected to be smaller [H since there is additional CKM suppression and 
color suppression. For the charged mode the corresponding tree diagram is external and not color 
suppressed. 




(c) (d) 



Figure 1: Feynman diagrams describing (a) B — B mixing; the decay B^ rj'K^ via (b) color- 
suppressed tree, (c, d) internal gluonic penguin. 

Additional higher-order amplitudes carrying different weak phases would lead to deviations, 
AS", between the measurements of the time-dependent CP violating parameter in these rare decay 
modes and in the charmonium decays. Theoretical bounds for these deviations have been 
calculated with an SU(3) analysis jHlIZl- Such bounds have been improved by recent measurements 
of B^ decays to a pair of neutral charmless light pseudoscalar mesons . From this and other recent 
experimental measurements, improved model-independent correlated bounds in the {S, C) plane 
for the decay B^ rj'K^ have been derived [H], with the conclusion that AS" is expected to be less 
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than 0.10 (with a theoretical uncertainty up to ~30% due to the assumptions in the calculation). 
Specific model calculations conclude that AS is even smaller; for instance a recent calculation 
jlOj finds AS* = 0.011 it 0.009 it 0.010, where the first error is due to theoretical uncertainties for 

— > rj'K^ and the second is for uncertainties in the charmonium-i^g system. A significantly larger 
value of AS" could arise from phases from non-SM amplitudes |3| . 

The CP-violating asymmetry in the decay — > rj'K^ has been measured previously by the 
BABAR jllj and Belle jl2j experiments. The measurement presented in this paper is an update 
of the previous BABAR measurement. In the present analysis the measurements of time-dependent 
CP violating parameters in B+ ij'K^ are used as a null control sample for the corresponding 
measurements in B^ —>■ r]'K^. 

2 The BaBAR Detector and Dataset 

The results presented in this paper are based on data collected in 1999-2004 with the BABAR 
detector ^31 at the PEP-II asymmetric e~^e~ collider Jl] located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 205 fb^^, corresponding to about 227 million BB pairs, was 
recorded at the T(4S') resonance ( "on-resonance" , center-of-mass energy ^/s = 10.58 GeV). 

The asymmetric beam configuration in the laboratory frame provides a boost of f3^ = 0.56 to 
the T{4:S). Charged particles are detected and their momenta measured by the combination of 
a silicon vertex tracker (SVT), consisting of five layers of double-sided detectors, and a 40-layer 
central drift chamber, both operating in the 1.5 T magnetic field of a solenoid. The tracking system 
covers 92% of the solid angle in the center-of-mass (CM) frame. 

Charged-particle identification (PID) is provided by the average energy loss (dE/dx) in the 
tracking devices and by an internally reflecting ring- imaging Cherenkov detector (DIRC) covering 
the central region. A K/ir separation of better than four standard deviations (a) is achieved for 
momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta in the B decay final states. 
Photons and electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The EMC 
provides good energy and angular resolutions for detection of photons in the range from 30 MeV to 
4 GeV. The energy and angular resolutions are 3% and 4 mrad, respectively, for a 1 GeV photon. 

The flux return for the solenoid is composed of multiple layers of iron and resistive plate cham- 
bers for the identiflcation of muons and long-lived neutral hadrons. 

3 Event Selection and Analysis Method 

Monte Carlo (MC) simulations of the signal decay modes, BB backgrounds, and detector response 
are used to establish the event selection criteria. 

We reconstruct B meson candidates by combining a Kg or a with an r]' meson. We select 
Kg — > TT+vr^ decays by requiring the invariant tt+vt" to be within a mass window of 12 MeV around 
the nominal K^ mass and requesting a flight length > 3 a. We select Kg tt^tt^ decays requiring 
the invariant vr^^vr'^ to be within a mass window of 30 MeV around the nominal Kg mass and a 
fit of the two vr'^ mesons to a common decay vertex. We reconstruct t]' mesons through the decays 
r]' — > p^j and r]' rj-K^iT^ with 7/ — > 77 or ^ tt'^tt'tt^. The photon energy must be greater 
than 50 (30) MeV for rj (t:^) candidates, and greater than 100 MeV in rj' — > p^j. We make the 
following requirements on the invariant mass (in MeV): 490 < m^^< 600 for r/^^, 120 < m^^ < 150 
for 7r° (100 < m^^ < 155 in Kg ^ 7r°7r0), 510 < m^^ < 1000 for p^, 520 < m^^^ < 570 for 773^, 
945 < m^f < 970 for t]' r^vr+vr^, and 930 < m^' < 980 for ??' p^j. 
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We make several particle identification requirements to ensure the identity of the signal pions. 
In charged B decays for the track we require an associated DIRC Cherenkov angle between 
—5a and +2 a from the expected value for a kaon. 

A B meson candidate is characterized kinematically by the energy-substituted mass m-ES = 
{^s + Po • pb^/Eq — and the energy difference AE = E*^ — \\fs-, where the subscripts and 
B refer to the initial T{AS) and to the B candidate, respectively, and the asterisk denotes the 
r(45) frame. We require |A^| < 0.2 GeV and 5.25 < toes < 5.29 GeV. 

To reject background in continuum e"^e~ — > qq events {q = u, d, s, c) , we make use of the angle 
6t between the thrust axis of the B candidate and that of the rest of the tracks and neutral clusters 
in the event, calculated in the center-of-mass frame. The distribution of cos 9t is sharply peaked 
near ±1 for combinations drawn from jet-like qq pairs and is nearly uniform for the isotropic B 
meson decays; we require |cos0t| < 0.9. The remaining continuum background dominates the 
samples and is modeled from sideband data for the maximum likelihood fits. 

We use Monte Carlo simulations of B^B^ and B~^B~ pair production and decay to look for BB 
backgrounds. Prom these studies we find evidence for a small BB background component for the 
channels with r/' — > p^j, and we have added a single BB component to the fit. 

Prom a BB pair we reconstruct a B^ decaying into the final state / = r]'Kg (Bcp)- We also 
reconstruct the vertex of the other B meson (i?tag) and identify its flavor. The time difference 
At = tcp — itagj where tcp and itag are the proper decay times of the CP and tagged B mesons, 
respectively, is obtained from the measured distance between the Bcp and Stag decay vertices and 
from the boost {(3^ = 0.56) of the e~^e~ beam system. The distribution of At is: 

e-|At|/T 

F{At) = — {1 =F Aw ± (1 - 2uj) [Sf sin(AmrfAt) - C/ cos(AmdAt)]} , (1) 

where the upper (lower) sign denotes a decay accompanied by a B^ (B^) tag, r is the mean B^ 
lifetime, Arud is the mixing frequency, and the mistag parameters uj and Au; are the average and 
difference, respectively, of the probabilities that a true B^ (B^) meson is tagged as B^ (B^). The 
tagging algorithm, based on six tagging categories, is an improved version of what was used in the 
previous BABAR publication Separate neural networks are trained to identify primary leptons, 
kaons, soft pions from D* decays, and high-momentum charged particles from B decays. Each 
event is assigned to one of the six mutually exclusive tagging categories based on the estimated 
mistag probability and on the source of tagging information. 

4 Maximum Likelihood Fit 

We use an unbinned, multivariate maximum-likelihood fit to extract signal yields and CP violating 
parameters. We indicate with j the species of event: signal, continuum background, or BB back- 
ground (ry' — > p^j). We use four discriminating variables: ?71es; AE, At, and a Pisher discriminant 
J^. The Pisher discriminant combines five variables: the angles with respect to the beam axis of the 
B momentum and B thrust axis in the T(4S') frame, the zeroth and second angular moments of 
the energy flow excluding the B candidate around the B thrust axis, and the tagging category. Por 
each species j and each tagging category c, we define a total probability density function (PDF) 
for event i as 

^ rjiruEs') ■ V,{AE') ■ V,{r) ■ V,{At\ aXi, c) , (2) 

where is the error on At for the event i. With rij defined to be the number of events of the 
species j and fj^c the fraction of events of species j for each category c, we write the extended 



10 



likelihood function for all events belonging to category c as 

= jY J H("'sig/sig,c'Psig,c + '^qgfm-.cPqq + ''^ B B fsig,c'P ^ ^ 
i 

where is the total number of input events in category c. In the last term of this formula we have 
assumed /sig,c = Jbb c- The total likelihood function for all categories is given as the product over 
the seven tagging categories (including a category for untagged events for yield determinations). 
The product is extended to additional sub-decays by multiplying the above product likelihoods for 
each decay. 

We maximize the likelihood function while varying a set of free parameters: 5, C, three back- 
ground J- PDF parameters, and, for each sub-decay, signal and background yields, Ai? and rriES 
background parameters, and six parameters representing the background At shape. 

5 Results 

The results of the fits to the neutral modes are shown in Table ^ We combine the five sub-decay 
modes shown in Table H The decay mode — > 'q'^^^^^K^{-K^TT^) has not been reconstructed 
because of its expected low signal yield. The results for the three charged modes, used as a 
crosscheck, are presented in Table |21 the values of S and C are consistent with zero, as expected. 
The first two columns in each table represent the primary modes that we reported on in our previous 
analysis. The last three columns of Table ^ and the last column of Table 121 are for decay modes 
that we have not reported on previously. The efficiency of the selection after all cuts is about 25% 
for the four primary decay modes and 15% for the four new decay channels. 

Table 1: Results from fits to sub-decay modes of — > r]'K^. 















Signal yield 


192 ± 15 


438 ± 27 


55 ±9 


50 ±9 


86 ±20 


BB yield 




93 ±33 






78 ±23 


S 


0.05 ±0.28 


0.41 ±0.19 


0.53 ±0.49 


-0.18 ±0.50 


-0.26 ±0.61 


c 


-0.12 ±0.18 


-0.29 ±0.13 


0.18 ±0.38 


-0.69 ±0.40 


0.19 ±0.44 


Combined fit 


Signal Yield 






819 ±38 






S 






0.27 ±0.14 






C 






-0.21 ±0.10 







In Fig. 12 we show projections onto rriEs and AE of a subset of the data of the primary decay 
modes for which the signal likelihood (computed without the plotted variable) exceeds a mode- 
dependent threshold that optimizes the sensitivity. In Fig. IHl we show the same projections onto 
TTiES and AE for the new decay modes considered in the present analysis. 

We show in Fig. |1] the At projections and asymmetry of all combined neutral modes for events 
selected as for Fig. |21and Fig. |21 In Fig. [S] we show the — 21n£ scans for S and C. The best fit 
value for S is 3.0 standard deviations from the BABAR value of sin 2/3 in charmonium decays [15j . 
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(3) 



Table 2: Results of fits to sub-decay modes rj'K^. 





rl 

'/jjTrTT-"- 




n' K+ 


Signal yield 


585 ± 26 


1322 ± 48 


221 ± 17 


BB yield 




774 ± 75 




S 


-0.13 ±0.13 


-0.05 ±0.10 


-0.31 ±0.23 


c 


0.00 ±0.10 


-0.10 ±0.08 


-0.06 ±0.16 


Combined fit 


Signal Yield 




2124 ± 57 




S 




-0.10 ±0.07 




C 




-0.05 ±0.06 






iTigg (GeV) AE (GeV) 

Figure 2: The B candidate mss and IS.E projections for B~^ — > r]'K~^ (a, b) and B^ r]'K^{c, d) 
in main decay modes . Points with errors represent the data, solid curves the full fit functions, and 
dashed curves the background functions; the shaded histogram represents the r]' K subset. 
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Figure 3: B candidate rriEs and AE projections for r]'K^{—^ vr^vr'^) (a, b), —>■ r]'^^^^-^^^Kg{^ 
vr+vr^) (c, d), and B^ r/^^g^-j^^K"*" (e, f). Points with errors represent the data, sohd curves 
the fuh fit functions, and dashed curves the background functions. The shaded histogram in (a, b) 
represents the rj'^^^^^^^K'^ subset. 



6 Systematic Uncertainties and Crosschecks 

The contributions to the systematic uncertainties in S and C are summarized in Table El We 
evaluate the uncertainties associated with the PDF shapes by variation of the parameters describing 
each discriminating variable. Systematic errors associated with signal parameters {At resolution 
function, tagging fractions, and dilutions) are determined by varying their values within their errors. 
Uncertainties due to Am^ and tb are obtained by varying these parameters by the uncertainty in 
their world average values ^16 . All changes are combined in quadrature obtaining an error of 0.01 
for both S and C. 

We vary the SVT alignment parameters in the signal Monte Carlo events by the size of mis- 
alignments found in the real data, and assign the resulting shift in the fit results as the systematic 
error. 

The systematic errors due to interference between the CKM-suppressed b ucd amplitude 
and the favored b — > cud for some tag-side B decays is found to be negligible for S and gives 
a contribution to the C uncertainty of about 0.012. The effect of BB background is estimated 
to be negligible, but we assign an uncertainty of 0.01 in S due to statistical limitations of this 
statement. An uncertainty of 0.02 is assigned to account for limitations of Monte Carlo statistics 
and modeling of the signal. We assign an uncertainty of 0.01 to account for the uncertainty in the 
position and size of the beam spot, determined from variation of these quantities in signal MC. 
The total systematic error is obtained by summing individual errors in quadrature. 

We have also performed a number of checks of our results. We divide the sample into two 
sub-samples — the previously published sample and the data collected in 2003/2004. We find 
consistency with our previous results and between the two sub-samples. When we fit with the 
value for C fixed to zero, we find S = 0.29 it 0.14. We produce samples of pseudo-experiments 



13 




Figure 4: Projections onto At for v'^s^ showing the data (points with errors), fit function 

(soHd hne), and background function (dashed hne), for (a) and (b) B^ tagged events and (c) 
the asymmetry between B^ and B^ tags. 
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s c 



Figure 5: — 21n£ scan for S (left) and C (right) parameters. The sohd blue line refers to all 
combined neutral sub-decays, the dotted line to the -q'^^^K^ sub-decay and the dashed line 

to the — > rjp^Kg sub-decay. The shaded band shows the BABAR value of sin 2/3 in charmonium 
decays [15] . 



Table 3: Estimates of systematic errors. 



Source of error 


a{S) 


aiC) 


PDF Shapes 


0.01 


0.01 


SVT alignment 


0.01 


0.01 


Tag-side interference 


0.00 


0.01 


BB Background 


0.01 


0.00 


MC statistics/modeling 


0.02 


0.02 


Beam spot 


0.01 


0.01 


Total 


0.03 


0.03 
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generated with events produced to match the PDF distributions. From these samples , we verify 
that the fit bias on S and C is neghgible and that there is a good agreement between expected and 
observed errors. 

7 Conclusion 

We have reconstructed 2124 it 57 r]'K~^ events and 819 it 38 v'^s events, about 

two-thirds of which have a flavor tag. We have used the latter sample to measure the time- 
dependent CP-violating parameters for B^ — > rj'K^. We find S = 0.27ib0.14(stat) ib0.03(syst) and 
C = -0.21 ± O.lO(stat) ± 0.03(syst). 

The measured value of S is 3.0 standard deviations from the the BABAR measurement of sin2/3, 
0.722 lb 0.040 lb 0.023, from B charmonium Kg decays |15j . The observed deviation AS is large 
compared with the expected theoretical uncertainty. 
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